function res = central_der2_onestep(u,N,x_max)
h = x_max/N;
u = fill_ghostcell_period(u,N,1);
res = (u(1:N+1)-2*u(2:N+2)+u(3:N+3))/(h*h); 